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g ; RESUMEN 

f-*) ' Se ha tratado de asociar a los rayos cosmicos de ultra alta energi'a (RCUAE) 

CSl ' a las fuentes mas variadas y poderosas que se conocen en el universe. Los 

! estallidos de rayos gamma (ERG) son candidates naturales. En este trabajo 

sostenemos que tales rayos cosmicos pueden ser acelerados por ondas elec- 
tromagneticas de alta amplitud (OEAA), las cuales son creadas al momento 
en que la aproximacion MHD del campo en el viento generado por el gener- 
ador central del ERG deja de ser valida. En tal situacion se considera como 
. generador central a una estrella extrafia con rotacion diferencial la cual ob- 

tiene la estrella al momento en que una estrella de neutrones se transforma en 
O ' una estrella extrana dentro de un sistema binario de rayos X de baja masa. 

^ , Las OEAA generadas de esta manera aceleran los iones ligeros a las mas al- 

I tas energias, E = qrjA^max, en donde la eficiencia rj ~ 10~^ da cuenta de 

todas las posibles perdidas energeticas. De manera alternativa, tambien con- 
sideramos la posibilidad de que una vez formadas, las OEAA se vuelquen 
^ ' entre si, dando paso asi a una turbulencia electromagnetica relativisticamente 

poderosa. Bajo esta suposicion, se estima un Ifmite inferior para la eficiencia, 
T] ~ 10~^-^. Debido a su edad, los sistemas binaries de rayos X de baja masa 
' pueden localizarse en regiones de baja densidad de materia interestelar, p. ej., 

. en ciimulos globulares o inclusive en el medio intergalactico en el caso de los 

I sistemas con un gran movimiento propio; en cuyo caso las perdidas energeticas 

sufridas por los rayos cosmicos en el medio circundante son evitadas y dichas 
particulas pueden experimentar el voltage maximo disponible, produciendo 
eventos que rebasarian inclusive las observaciones existentes. 
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^ ■ ABSTRACT 



Ultra-high-energy cosmic rays (UHECRs) have been tried to be related to 
the most varied and powerful sources known in the universe. Gamma-ray 
bursts (GRBs) are natural candidates. Here, we argue that cosmic rays can 
be accelerated by large amplitude electromagnetic waves (LAEMWs) when 
the MHD approximation of the field in the wind generated by the GRB's 
magnetized central engine breaks down. The central engine considered here 
is a strange star born with differential rotation from the accretion induced 
conversion of a neutron star into a strange star in a low-mass X-ray binary 
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system. The LAEMWs generated this way accelerate hght ions to the high- 
est energies E = qrjA^rnax with an cfRcicncy t] ^ 10~^ that accounts for all 
plausible energy losses. Alternatively, we also consider the possibility that, 
once formed, the LAEMWs are unstable to creation of a relativistically strong 
electromagnetic turbulence due to an overturn instability. Under this assump- 
tion, a lower limit to the efficiency of acceleration is estimated to be about 
r] ~ 10^^-^. Due to their age, low mass X-ray binary systems can be located 
in regions of low interstellar medium density as, e.g., globular clusters or even 
intergalactic medium in case of high proper motion systems, and cosmic ray 
energy losses due to proton collisions with photons at the decelerating region 
are avoided, thus opening the possibility for particles to exploit the full voltage 
available, well beyond that currently observed. 

Key Words: Cosmic Rays — Gamma-ray sources — gamma-ray bursts — 
Neutron stars 

1. INTRODUCTION 

A long standing puzzle in astrophysics remains the origin of the highest 
energy cosmic ray events with energies in excess of 10-^° eV (Nagano & Watson 
2000). The ultra- high-energy cosmic rays (UHECRs) observed so far follow 
a well known energy spectrum N{E) ~ E~'^-'^ beyond the feature known as 
the 'ankle' located at li]^^-^ cV (Axford 1994: Bird et al. 1994). The origin 
of the particles with energies greater than ~ 10^^ eV is usually attributed to 
extragalactic sources (Axford 1994). The distance an UHECR can travel is 
limited by the energy losses due to photo-pion production in collisions with 
cosmic microwave background photons. This sets a radius, Rgzk ^ 50 Mpc, 
beyond which the energy losses are copious and where the energy of a cosmic 
ray can be lowered even by a few orders of magnitude. This mechanism would 
generate a signature in the spectrum at ~ 5 x 10^^ eV, known as the GZK 
cutoff (Greisen 1996; Zatsepin & Kuzmin 1966). 

However, in recent years serious constraints on its existence prompted the 
attention to models containing new physics as the possible solution for the 
UHECR enigma (see Bcrczinsky 2008 for a review of the models). 

A plethora of astrophysical objects have been considered as possible sources 
of UHECRs. Among them gamma ray bursts (GRBs) have received special 
attention (Vietri 1995; Milgrom & Usov 1995; Waxman 1995; Chen, et 
al. 2002; Vietri et al. 2003). In fact, the energetics associated to local 
GRBs is found to enable a GRB-UHECR association (Wick et al. 2004). 
Wick et al. (2004) considered that a good estimate for the local emis- 
sivity el^HECR needed to power cosmic rays with energies > 10^° should be 
euHEcn « U/t^^ « lO'^'^C/ai ergs Mpc^^ yr-\ where U = lO^ij/ai ergs cm'^ 
is the observed energy density in > 10^° eV cosmic rays {U21 « 0.5 for the 
HiRes experiment and U21 w 2 for AG AS A) and the timescale of a proton with 
energy ^ 10^° cV for photo-pion production is t^^. For a local distribution 
of GRBs (z <IC 1), which is inside the GZK cutoff, pre-BeppoSAX estimates 
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of the luminosity density (Vietri 1995; Waxman 1995) argue for a value of 
t'cRB^ ~ 10^^644 ergs Mpc~"^ yr""'^, where £44 « few. In a more recent analysis, 
Vietri et al. (2003) have arrived at a value of £44 « 1.1. Even considering 
the uncertainties in these estimates, the similarity of the two luminosity den- 
sities e'-^HECR and k'^cRs^ make it reasonable to consider that UHECRs may 
be powered by the local GRBs. 

The ~ 10^° - 10^2 ergs (Frail 2004) required to fit the estimated GRBs' 
energetics, as well as the short observed GRB durations, naturally point to- 
ward gravitational energy associated to compact objects as the original en- 
ergy source. One broad class of central engine (CE) models consider the 
merger of two neutron stars (Paczyhski 1986), or a black-hole and a neutron 
star (Paczynski 1991), in a binary system or an isolated "failed" supernova 
(Wooslcy 1993) in which the released gravitational energy is temporarily 
stored in thermal energy and part of it, through neutrinos, is inefficiently con- 
verted into relativistic outflow and successive shocks required to fit the GRBs' 
light curves. Most of these models eventually involve a black-hole which also 
swallows a very large part of the energy. And in addition to that, the recent 
analysis of (Soderberg et al. 2006) shows that the supernova-GRB connection 
is not as strong as previously tought, inducing (Gal- Yam et al. 2006) to sug- 
gest a new explosive process to explain the lack of GRB-associated supernova 
detections. Another class of CE models takes advantage of the possibility of 
storing a significant part of the released gravitational energy into rotational ki- 
netic energy of the newly formed compact object (Usov 1992), hence limiting 
the tremendous energy drain by thermal neutrinos and avoiding the black-hole 
gluttony. This kinetic energy can then be extracted on time scales of seconds 
using the torque exerted by the strong magnetic field of the compact star. In 
the "DROCO" (Differentially Rotating Compact Object) scenario (Kluzniak 
& Ruderman 1998) differential rotation of the compact object acts as a pow- 
erful dynamo and results in magnetic field strengths of the order of 10^^ G or 
higher where subsequent buoyancy forces induce emergence of flux ropes at 
the stellar surface and produce successive sub-bursts. 

2. DIFERENTIALLY ROTATING COMPACT OBEJCTS 

Here, we explore the possibility that UHECRs are accelerated in the 
Poynting-dominated outflow of a GRB when the central engine is a new born 
strange star with a rotation period ~ lO'* s~^ and a magnetic field of the 
order of 10^^ G. Strange stars (Alcock et al. 1986; Haensel et al. 1986) are 
self-bound compact object made entirely of deconfined quark matter, com- 
prising the three flavors of light quarks u, d and s, which is assumed to be 
the true ground state of hadronic matter. Within this "strange matter hy- 
pothesis" (Boclmcr 1971; Wittcn 1984) once a seed of quark matter appears 
within a compact star it will rapidly convert the whole star into a strange star 
in a time scale of the order of milliseconds up to a few seconds (Lugones et 
al. 1994; Harko et al. 2004). If this can happen in core collapse supcrnovae 
then all objects we call "neutron stars" may actually be strange stars; there 
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is however some evidence that this is not the case (see, e.g., (Weber 2005) 
for a discussion). We will consequently consider a more plausible scenario 
assuming that the critical density for appearance of quark matter is large 
enough that only the most massive neutron stars may convert into strange 
stars. Such a conversion is thus expected to occur in neutron stars which 
have accreted a significant fraction of a solar mass in a low-mass X-ray bi- 
nary system (LXMB). The energetics, comparable to a typical core-collapse 
supernova, and the expectable frequency of such events were considered by 
Cheng & Dai (1996) who proposed them as candidates for GRBs. Given the 
accretion spin-up of the neutron star and its subsequent shrinking during its 
conversion to a strange star, the new-born strange star is plausibly rotating at 
sub-millisecond period while the differential rotation due to the very different 
density profile of a strange star compared to a neutron star can be expected 
to results in a DROCO as described in (Dai & Lu 1998), i.e., with a magnetic 
field of the order of 10^^ G. An attractive aspect of the strange star scenario 
is that, most baryonic matter being converted into quark matter, it avoids the 
problem of baryon contamination of the outflow (Haensel ct al. 1991). Nev- 
ertheless, the outer crust of the neutron star, with a mass of about 10~^ Mq, 
is most probably not converted into strange matter and rather expelled by 
the neutrino flux and the highly super-Eddington thermal luminosity of the 
bare strange star (Usov 1998, 2001; Page & Usov 2002), as well as by the 
flux ropes that emerge at the stellar surface (Kluzniak & Ruderman 1998), 
providing the protons which can be accelerated at energies well above 10^'^ as 
we will describe later. 



3. THE FLUX 



The particle flux which one can expect from such event can be estimated 
by analogy with standard pulsar physics. In the magnetosphere of the CE 
a number of processes take place that enhance the multiplicity of pairs com- 
pared to the proton current (Usov 1992), besides the pairs produced by the 
strange star surface. For an estimate we consider the pulsar's proton current 
as determined by the Goldreich-Julian current (Goldreich & Julian 1969) 

N, = f Ngj^ f^^^ -10'^ ^If B,, s-i (1) 
ec 

where Q.^ = Q/IO* s~^ and Bn = B/IO^"^ G, and have introduced an en- 
hancement factor / which could be larger than 1. We can write the particle 
flux at the Earth, for energies > 10^" eV, as 

/(> lO^OeV) ~ N,At auoHB naT2o (2) 

where At is the duration of acceleration activity of the source, i^crb the local 
GRB rate, a fraction a of them being of the type we consider, Ua = 0.02 
Mpc~^ the local density of galaxies and T20 the lifetime of protons above 
10^" eV. For energies in excess of 10^" eV, correspond propagation times 
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T20 < 3 X 10^ yr, however, the vast majority of the sources within the region 
Hmited by Rqzk lie at '--^ 40 — 50 Mpcs from us, thus it is appropiate to 
use T20 ~ 3 X 10^ yr (Aharonian & Cronin 1994). So far the true birthrate of 
local GRBs has not been established. The enhancement factor which accounts 
for the coUimation of the luminosity has been calculated recently by a fitting 
procedure to both the observed differential peak flux and redshift distributions 
at the same time (Firmani et al. 2004), yielding a local GRB rate Vgrb ~ 
5 X 10~^ yr~^. A reasonable estimate (see below) for a is of the order of 10%. 
This gives us, considering a production life-time At ~ 10 s, 

/(> lO^^eV) ~ 2 X 10-2VSl7^^lA^lo cm-2 s"^ sr'^ (3) 

which is safely comparable, and even above if / > 1, to the observed flux 
/(> lO^OeV) ~ 2.6to;6 X 10"^° cm'^ s'^ sr"! (Nagano & Watson 2000). 
Notice that the energy in cosmic rays of 10^° eV, 7Vp(At)10^°eV ^ lO^'^ ergs, 
represents about 1% of the total initial kinetic energy of the star. Since not 
all the kinetic energy lost by the star will necessarily be detectable, in the 
GRB, this energy may be comparable, or even larger, than the GRB energy. 

Acceleration of cosmic rays in GRBs has been considered, as in many other 
astrophysical objects, to be due to the Fermi acceleration process (Fermi 
1949), this mechanism taking place either in the plasma outflow (Waxman 
1995) or in the external shock of the outflow when the latter hits the inter- 
stellar medium at a large distance from the CE (Vietri 1995). Nonetheless, 
it has been shown by Gallant and Achterberg (Gallant & Achterberg 1999) 
that the energy gain by this mechanism in external shocks when starting from 
nonrelativistic energies is much lower than predicted previously, so in a more 
realistic situation the Fermi acceleration process could hardly reproduce the 
highest energy events so far observed. Chen, et al. (2002) outlined a different 
scenario in which Alfven shocks traveling through the plasma can accelerate 
light nuclei in a similar fashion a terrestrial linear acceleration does, i.e., the 
acceleration is stochastic and particles are accelerated by the interplay of the 
accelerating and decelerating phases of the plasma wakefields generated by 
the Alfven shocks. In this stochastic process, the energy losses associated 
with curvature and synchrotron radiation and with ultra-high-cnergy proton- 
plasma proton collisions are minimal, so the system can be regarded as one 
in which the acceleration is maximum. 

Instead of considering Alfven shocks, we consider here that large-amplitude 
electromagnetic waves (LAEMWs) are responsible not only of accelerating 
particles to produce a GRB event (Usov 1994) but also responsible for the 
acceleration of plasma protons to the highest energies. Poynting flux domi- 
nated models of GRBs have been explored recently (Zhang & Meszaros 2004) 
as an alternative to the fireball model of Recs & Meszaros (1992) because the 
efficiency of radiation is rather high and the peak photon energy distribution 
Ep could be well reproduced, among other observational facts. LAEMWs are 
naturally generated when the MHD approximation of the magnetized out- 
flow breaks down. This happens at a radius, Rmhd, from the CE where 
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the plasma density n±, which decreases hke ~ (Rhoads 1999), drops 

below the density needed for applicability of the MHD approximation, the 
Goldreich-Julian density Ngj, which decreases like ~ r"^ (Goldreich & Ju- 
lian 1969). The LAEMWs generated this way have a frequency f2 ~ 10^ s~^, 
change polarity on a length scale of 27r(c/f2) ~ 2 x lO'' cm, and are very 
strong which can be seen if we consider the adimensional Lorentz-invariant 
parameter u = eB/mcfl ~ 10^ which is the ratio of the cyclotron frequency 
in the wave field w = eB/mc to the wave frequency (Gunn & Ostriker 1971). 
Such a high ly value can be interpreted by saying that the particles move in a 
nearly static field. Particles surf-ride (Arons 2003) the LAEMWs and suffer a 
ponderomotive force in the radial direction due to the inductive electric field 
orthogonal to the magnetic field fpond = {<l/c){Sv x SB) = mcflv. Because the 
particles move relativistically, the radius of curvature associated to an UHE 
proton fc = 2i?;r^ ^ r, where Ri = mc^-f/eB is the Larmor radius and is 
the wind's Lorentz factor (Arons 2003), and hence excessive synchrotron and 
adiabatic losses are negligible. In principle, following Ruderman & Sutherland 
(1975), the maximum differences in the voltage of the wind outflow at large 
distances from the central engine are 

2c2 V 10'* s- 

« ^ V. (4) 



2 X 1017 G / V 10^ cm 



so that we have an extractable energy well beyond the currently observed 
data. The work produced by the ponderomotive force can be approximated 
by 

max (5) 

where r] is the efficiency factor. From analysis of pulsar wind nebulae it is seen 
that about 10% of the total voltage is transferred to the ions (Arons 2003). 

In the situation prevailing in a GRB outflow, an UHE proton can loose 
some of its energy due to collisions with ambient protons. This can be esti- 
mated in an easy way (Chen, et al. 2002) by integrating the collision proba- 
bility, UppHpir) , from Rmhdi the transition radius, to Rmhd~\~ Rmfp where 
Rmfp is the mean-free path radius and the collision cross section Upp 30 
mb. If we want the system be collision-free, the threshold condition 



ctppNgjRmhd 



< 1 (6) 



should be satisfied, where we have taken np{r) = Nqj for the plasma density 
at the transition region. It is easy to see that for appropriate values this 
condition is always satisfied. 

However, as an alternative to this high efliciency acceleration mechanism, 
we have also considered the possibility that the interaction of LAEMWs with 
the plasma can lead to instabilities. Lyutikov & Blackman (2001) argued that 
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the interaction of LAEMWs with the plasma can lead to violent instabilities 

where the MHD approximation breaks down because an oveHurn instability 
of LAEMWs, creating a broad spectrum of random electromagnetic fields. 
The mixing of various parts of the wave destroys its oscillatory structure 
thus yielding a relativistically strong electromagnetic turbulence. When the 
correlation length, Ic, of the turbulent electric field is much shorter than the 
radiative length, 1^ Ir, there will be an effective root-mean-squared deficit 
in the coherence of the ion velocity and the accelerating electric field. If 
the turbulent energy is concentrated near the smallest scale ~ skin depth, a 
lower limit to the efficiency of acceleration, rj, is obtained, which in case of 
protons can be expressed in terms of the classical proton radius rp and the 
wave frequency fl, rj = (2rpn/c)i/^ » IQ-^-^ We can see that even if this 
lower factor is taken into account, we have an extractable voltage sufficient to 
account for the highest energy events. 

Whether the turbulence is considered or not, the principal energy losses 
for a high energy proton once the acceleration process is almost completed 
are due to collisions with the intense thermal photon field left after adiabatic 
expansion of the relativistic outflow. Each collision undermines ~ 10% of the 
UHE proton energy (Aharonian & Cronin 1994) although the characteristic 
time for this energy loss mechanism (Waxman 1995) tpj ^ 10 /n^apjC is 
larger than the time of expansion of the outflow even when this is entering the 
decelerating regime. The collision cross section in this case is ap.y « 10~^^ cm^. 
The photon number density n-y was calculated by using the photon luminosity 
in the decelerating region, = A-Kr^c^n-yC-y, where e-y is the observed photon 
energy, i.e., ~ 1 MeV. 

Due to their age, LMXB systems are not associated with star formation re- 
gions and an interesting possibility arises if we consider the fraction of LXMB 
systems related to globular clusters, ^ 10% (Liu et al. 2001), in wich the 
CE is surrounded by a very tenous ambient medium. We would thus expect 
in this case that an equal fraction of the GRBs discussed here would not pro- 
duce an afterglow, or at least a very faint event (Castro-Tirado et al. 2005). 
This is another feature which avoids energy losses of UHECRs due to the 
photon field of the afterglow. However, we cannot exclude that the impact 
of the burst on the strange star's companion will produce a detectable signal 
(Ramirez-Ruiz 2004), but with negligible effect on the bulk of UHECRs since 
it will be of small angular size, restricted to the direction of the companion, 
and, in addition, the zone where the acceleration takes place is far beyond the 
location of the companion. 

Details of the expected UHECR spectrum are impossible to predict with 
confidence. In case of pure electromagnetic spindown of the central engine 
(Arons 2003) an _E^^ dependence is expected while if acceleration can be 
regarded as a quasi-stochastic process a spectrum ~ E^'^ would be more 
appropriate (Chen, et al. 2002). Propagation effects are usually argued to 
cause a steepening of the spectrum, and to that respect it has been argued 
that the consistancy between the data obtained by the HiresI-II, Yakutsk, 
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and Auger experiments proves the existence of a GZK cutoff in oposition 
to the results obtained by the Akcno-AGASA experiment, (Bcrczinsky 2008; 
Kampert 2008), but the lack of robust data makes impossible to reach a final 
conclusion. The model as presented here would agree with an extension of 
events beyond the highest energy values observed, > 10^'^ eV. 

Finally we mention that the above outlined scenario can be applied equally 
well to the accretion induced collapse of a white dwarf into a neutron star, 
considering a cataclismic variable as an initial system instead of a low mass 
X-ray binary. However, a young strange star with a bare quark surface will 
emit, for a possibly long period of time, super-Eddington luminosities through 
pair production at its surface (Usov 1998, 2001; Page & Usov 2002) and later 
through electron-electron bremsstrahlung from its electrosphere ( Jaikumar et 
al. 2004), with unique spectral characteristics which would allow to identify 
it unambiguously. On the contrary, a neutron star born through accetion 
induced collapse of a white dwarf is difficult, if not impossible, to be differ- 
enciated from its sibblings born through the core collapse of a massive star. 
Identification of strange stars in binary systems (see, e.g., for some possibilities 
(Page & Gumming 2005)) would give support to our scenario. 

The author thanks V. V. Usov for comments on the manuscript. This 
work was supported by grants from the Mexican Conacyt (# 36632-E) and 
from UNAM-DGAPA (# INI 12502). 
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